Inhomogeneous Loop Quantum Cosmology: Hybrid Quantization of the Gowdy Model 
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The Gowdy cosmologies provide a suitable arena to further develop Loop Quantum Cosmology, 
allowing the presence of inhomogeneities. For the particular case of Gowdy spacetimes with the 
spatial topology of a three-torus and a content of linearly polarized gravitational waves, we detail 
a hybrid quantum theory in which we combine a loop quantization of the degrees of freedom that 
parametrize the subfamily of homogeneous solutions, which represent Bianchi I spacetimes, and a 
Fock quantization of the inhomogeneities. Two different theories are constructed and compared, 
corresponding to two different schemes for the quantization of the Bianchi I model within the 
improved dynamics formalism of Loop Quantum Cosmology. One of these schemes has been recently 
put forward by Ashtekar and Wilson-Ewing. We address several issues including the quantum 
resolution of the cosmological singularity, the structure of the superselection sectors in the quantum 
system, or the construction of the Hilbert space of physical states. 

PACS numbers: 04.60. Pp,04.60.Kz,98.80.Qc 



I. INTRODUCTION 

In the absence of a full theory of canonical quantum 
gravity, it is most instructive to study symmetry-reduced 
systems in order to test and develop mathematical tech- 
niques that can help in achieving the quantization of gen- 
eral relativity, as well as to progress in the understanding 
of the physical phenomena emerging in quantum gravity. 
The symmetry reduction makes these systems more man- 
ageable for a complete quantization and facilitates the 
extraction of physical predictions. In particular, mim- 
icking the techniques applied in Loop Quantum Grav- 
ity (LQG) [lj, a quantization program for homogeneous 
models, known as Loop Quantum Cosmology (LQ C) Q, 
has been developed in recent years (see e.g. |3l-[l3|L The 
so-called polymeric quantization applied to these sys- 
tems has interesting physical consequences. Remarkably, 
the analog of the classical cosmological singularity is not 
present in the quantum theory and is therefore resolved. 

In view of the success of LQC, it is compelling to ex- 
tend the application of its techniques to the quantization 
of inhomogeneous cosmological systems. Several reasons 
motivate this extension. On the one hand, this kind of 
models preserves one of the most characteristic features 
of the full theory, namely the field-like nature of the de- 
grees of freedom. Thus, one would expect that their 
quantization gives trustable insights about some of the 
open questions in LQG which are related to the pres- 
ence of an infinite number of degrees of freedom. On the 
other hand, it is essential to investigate the role played 
by inhomogeneities in the quantum theory in order to 
understand the physical laws that explain the origin and 
evolution of the Universe, and reach in this way realis- 
tic predictions that can be confronted with cosmological 
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observations. Furthermore, the analysis of the inhomo- 
geneities would allow us to check the robustness of the re- 
sults obtained in homogeneous LQC, in particular those 
concerning the quantum resolution of classical cosmolog- 
ical singularities. 

With this aim, we have carried out a complete quan- 
tization of the linearly polarized Gowdy T 3 model (14J . 
We have chosen this model because it is the simplest 
cosmological system in vacuo which contains inhomo- 
geneities. Indeed, not only its classical solutions are 
well known [la . Ilq |. but also their quantization (by 
more conventional methods) has deserved much atten- 
tion since the 70s [17| • Actually, a complete Fock quan- 
tization of the deparametrized system has been achieved 
[18| , which has been shown to be essentially unique [19J . 
Unlike these previous analyses, which involve standard 
non-polymeric techniques, we will discuss here a hy- 
brid quantization that combines the polymeric proce- 
dures of LQC — applied to the homogeneous sector of 
the system, namely the set of degrees of freedom that 
describe the subfamily of homogeneous solutions of the 
Gowdy model — with the Fock quantization of the inho- 
mogeneities. Since the full polymeric quantization of the 
model is an extremely complicated task, we have chosen a 
more conservative approach, which explores the effects of 
the quantum discrete geometry underlying LQC only on 
the homogeneous gravitational sector. A natural treat- 
ment for the inhomogeneities is then a Fock quantiza- 
tion. This approach assumes that there exists a regime 
in which the most relevant phenomena emerging from 
quantum geometry are those affecting the homogeneous 
subsystem, whereas such effects are small and can be ig- 
nored for the inhomogeneities, even if the latter may still 
present a standard (Fock) quantum behavior. Although 
our quantization approach is not completely derived from 
LQG, it retains interesting features associated with the 
discrete and polymeric nature of the geometry, declaring 
a certain type of perturbative hierarchy on their relevance 
for different subsectors of the cosmological system. 



As we have commented, the homogeneous sector of 
the Gowdy model coincides with the phase space of the 
Bianchi I model. It turns out that there exist two dif- 
ferent quantizations of this spacetimc in the literature of 
LQC, developed in Refs. [ll| and [12j], which correspond 
to two different schemes of the so-called improved dy- 
namics prescription and that, for brevity, we will denote 
as schemes A and B, respectively. In Refs. [2y, |2l[, we 
presented our hybrid quantization making use of scheme 
A in the representation of the ho mog eneous sector, i.e., 
applying the quantization of Ref. [ll| (see also Ref. (9||). 
The main aim of this paper is to revisit and develop 
our hybrid quantization by incorporating the alternative 
scheme B, put forward by Ashtekar and Wilson- Ewing, in 
the construction of the Bianchi I representation for the 
homogeneous sector. We will review the quantization 
of Ref. [12| , improve and complete some aspects of the 
quantization procedure, extend our hybrid framework to 
adapt it to the alternative scheme and, finally, compare 
the results obtained with the two hybrid approaches. 

With this objective, the constraints present in the 
model will be represented by well-defined operators on 
a (hybrid) kincmatical Hilbert space. These constraints 
are a global diffeomorphism constraint, which generates 
translations in the circle, and a global Hamiltonian con- 
straint which couples the homogeneous and inhomoge- 
neous sectors in a complicated form. As it is typical in 
the context of LQC, in both schemes the Hamiltonian 
constraint operator superselects the kinematical Hilbert 
space — which as a whole is non-separable — in separa- 
ble sectors. In particular we will be able to determine 
the structure of these superselection sectors for scheme 
B, which are highly non-trivial and had not been identi- 
fied previously (see Ref. [12|]). Furthermore, our Hamil- 
tonian constraint leads in fact to a difference equation 
in an internal, strictly positive parameter with support 
on semilattices of points with constant separation. In 
scheme A, this parameter is the classical global time used 
in the Fock quantization to deparametrize the system 
(l8| . In Refs. [2fJ, [21| we showed that the solutions to 
the Hamiltonian constraint are determined by the data 
on the initial section, namely the section where that pa- 
rameter takes its lowest allowed value. We will see that, 
in scheme B, the analogous discrete variable is (up to a 
constant factor) the volume of the Bianchi I spacetime 
associated with the Gowdy universe, and we will analyze 
the resulting structure of the solutions. 

Actually, one of the motivations of our work is the ex- 
pectation that the quantum field theory should remain 
valid on the loop quantized Bianchi I background after 
imposing the quantum constraints, therefore validating 
the standard Fock description for the inhomogeneities in 
the approximation adopted here. In the hybrid approach 
for scheme A, we demonstrated in Refs. [2(J, UM that 
the physical Hilbert space has the expected tensor prod- 
uct structure: the physical Hilbert space of the Bianchi I 
model times a Fock space for the inhomogeneous sector, 
which is equivalent to that obtained in Ref. [18]. Hence, 



we indeed recover the standard quantum field theory for 
the inhomogeneities. We will discuss whether this contin- 
ues to be the case in the new hybrid approach, obtained 
with the alternative quantization scheme for Bianchi I. 

Another important motivation, as we have commented, 
is the study of the fate of the cosmological singularities 
in the quantum theory. We will see that the polymeric 
quantization of the homogeneous sector, irrespective of 
the adopted improved dynamics scheme, is enough to 
cure the singularity. This fact endorses the previous re- 
sults of singularity resolution obtained in homogeneous 
LQC. 

The rest of the paper is organized as follows. In Sec. |TT] 
we summarize the basic features of the classical Gowdy 
T 3 model. In Sec. IIIII we construct the kinematical 
Hilbert space, and represent the elementary operators 
in both schemes A and B. In addition, we represent the 
diffeomorphism constraint in the quantum theory. We 
construct the Hamiltonian constraint operator in Sec. 
IIVI for the two studied schemes. In Sec. [V] we review 
the structure of the physical Hilbert space for scheme A. 
The physical sector for case B is investigated in Sec. IVII 
In Sec. IVII1 we discuss the main results of this work, 
and conclude with some additional remarks. Finally, we 
include four appendices, which contain some supplemen- 
tary technical details. 



II. THE REDUCED MODEL 

Gowdy models describe globally hyperbolic spacetimes 
in vacuo which possess two spacelike commuting Killing 
vector fields and whose spatial sections are compact [lj] . 
The simplest case is the Gowdy T 3 model (i.e., the spa- 
tial topology is that of a three-torus). We will focus 
our discussion on the subsystem with linearly polarized 
gravitational waves. This model satisfies the additional 
restriction that the Killing fields be hypersurface orthog- 
onal. Hence, they can be chosen mutually orthogonal ev- 
erywhere. Let d a and ds be these two axial Killing fields. 
We may then describe the spacetime (globally) with co- 
ordinates adapted to these symmetries: {i, 6*, cr, <5} , with 
0,cr,(5 £ S . The metric components depend only on t 
and 6. Since they are periodic in 6, they can be expanded 
in a Fourier series. 

We use the symmetries of the system to fix completely 
the gauge freedom associated with the diffeomorphism 
constraints in a and 5, the directions defined by the 
Killing fields, as explained in Ref. [22|]. We further fix 
all the gauge freedom associated with the inhomogeneous 
(non-zero) modes of the ^-momentum constraint and of 
the Hamiltonian constraint. The details of this gauge 
fixing can be found in Ref. [2l[ . At the end of the day, 
two global constraints remain in the model: a generator 
of translations in the circle, C'e , and a Hamiltonian con- 
straint, Cq. Except for these constraints, all the gauge 
freedom is fixed. The reduced phase space contains three 
pairs of "point-particle" degrees of freedom, which pa- 



rameterize the sector of (spatially) homogenous space- 
times and is therefore called the homogeneous sector, and 
all the non-zero modes of a field and its conjugate mo- 
mentum, which form the so-called inhomogeneous sector. 
In order to quantize the homogeneous sector with 
polymeric techniques according to LQC, we need to de- 
scribe it in terms of Ashtekar variables. The subfam- 
ily of homogeneous solutions of the Gowdy model, en- 
coded in this sector as we have mentioned, corresponds 
to empty Bianchi I spacetimes with three-torus topol- 
ogy. In a diagonal gauge, the non-trivial components 
of the SU(2) gravitational connection and of the densi- 
tized triad are Cj/(27r) and pi/(A-K 2 ), respectively, with 
{ci,pj} = 8-irGj6ij and i,j = 9, a, 6 (see e.g. [10i - fia| ) . 
Here G is the Newton constant and 7 the Immirzi pa- 
rameter. For the inhomogeneous sector we will carry out 
the Fock quantization presented in Ref. JJ| ■ With this 
aim we adopt as basic variables the creation and anni- 
hilation variables which would be naturally associated 
with a free massless scalar field, (a m ,aj„), where m can 
be any non-zero integer. In Appendix lAl we provide the 
form of the metric in terms of the chosen variables. The 
constraints are given by [21| 
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The first line of Cq is just the standard form of the den- 
sitized Hamiltonian constraint of the Bianchi I model in 
Ashtekar variables (see e.g. [10|). The inhomogeneities 
are present in the term Hq , which is the Hamiltonian cor- 
responding to a free massless scalar field, and in the term 
H^ nt , which represents an interaction term quadratic in 
the field. Note that the inhomogeneities are coupled 
with the homogeneous sector in a complicated way, what 
makes the hybrid quantization of this system far from 
straightforward. 



III. THE KINEMATICAL HILBERT SPACE 

The kinematical Hilbert space of this reduced Gowdy 
model is the tensor product of the kinematical Hilbert 
spaces of the two sectors. 



A. The homogeneous sector 

In order to construct the kinematical Hilbert space of 
the homogeneous sector, we consider the two different 
polymeric quantizations of the Bianchi I model developed 
in Refs. [ll|, Qjj]. The main distinction between them lies 
in the way in which the quantization prescription known 
as improved dynamics (introduced in Ref. [4j with suc- 
cessful results in isotropic cosmologies) is adapted to the 
anisotropic case. Although the quantization of Ref. [ll| 
is well defined in the considered model with compact spa- 
tial slices, it suffers from some drawbacks in non-compact 
situations [1(|. This fact motivated the consideration of 
an alternative prescription, whose corresponding quan- 
tization has been recently developed by Ashtekar and 
Wilson- Ewing in Ref. [T2|. In the following, we will 
call scheme A the prescription adopted in Ref. [ll|, and 
scheme B that of Ref. [I4I . 

Let us summarize the main characteristics of both 
quantizations. The elementary configuration variables 
are holonomies of connections computed along edges of 
oriented coordinate length 27r/Zi in the direction i, where 
fii is any real number. On the other hand, the elemen- 
tary momentum variables are triad fluxes through rect- 
angles orthogonal to those directions. The configuration 
algebra is the algebra of almost periodic functions gener- 
ated by the matrix elements of the holonomies, namely 
A/)i- (cj) = exp(i/ZjCj/2) (no Einstein summation conven- 



tion is adopted). We call Cylg 



iCyl s j the corre- 



sponding vector space, and employ the ket notation |//$ 
to denote the states jY^Cj) in the triad representation. 
The kinematical Hilbert space, "HKm = ^i'HKin,*) is then 
the completion of the space Cyl s with respect to the dis- 
crete inner product (/Xj|/4) = Sum'. f° r each direction. 
Hence, the states \fo), which are eigenstates of the op- 
erator pi associated with fluxes, provide an orthonormal 
basis of 'Hkum- In turn, the operators A/O associated 
with holonomies produce a shift equal to u' in the label 
Hi of this basis of states (see Refs. [3,[9|;[L1]])- We assume 
that the action of any operator defined on Cyl s is the 
identity when acting on the kinematical Hilbert space of 
the inhomogeneous sector, which will be introduced be- 
low. 



1. Improved dynamics: schemes A and B 

The so-called improved dynamics prescription is based 
on the assertion that, because of the existence of a min- 
imum gap A in the spectrum of the physical area oper- 
ator in LQG [J, [23[ , fiducial surfaces cannot be as small 
as one wants, but there is a minimum fiducial surface, 
Smm, whose corresponding physical area is equal to A. 
As a consequence, the coordinate length of the holonomy 
along each edge turns out to exhibit a minimum non- 
zero value 2-KfJLi. Prescription A is derived assuming that 
S'min, lying e.g. in the j — fc-plane, is a fiducial square 
with minimum fiducial side determined by 27r/ii, whereas 



according to prescription B it must be instead a rectan- 
gle with minimum fiducial sides given by 2irftj and 2-Kp^. 
Here, and in what follows, whenever the three indices i, 
j, and k appear together, we assume e^ 7^ 0. 

As a result one gets two different expressions for pi, 
one for each prescription, given by 
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3. Quantum representation in case B 

Similarly, one can introduce new parameters \ for the 
three fiducial directions such that the action of the op- 
erator Npn has only a non-trivial effect on the label Xi, 
whereas it does not change the other two labels of the 
state, Xj and A& [l2j. Remarkably, this effect is not 
a constant shift anymore but depends on the values of 
those two labels. Nonetheless, as we have commented, 
there exists a variable v, given in terms of the A^'s by 
v = 2AeA CT A,5, such that all the operators M^ produce a 
constant shift on it (for fixed orientation of the triad com- 
ponents). Therefore, it is convenient to work e.g. with 
the relabeled states \v, X a , Xs). The two A's are variables 
which measure the degree of anisotropy. The representa- 
tion of the basic operators is determined by [l2| 



Although these arguments can be considered heuris- 
tic inasmuch as the true relation between full LQG and 
LQC is still to be understood, there is a feature of pre- 
scription B (deduced in Ref. [12] following a procedure 
of this type) which distinguishes it and has not been re- 
alized until now in the literature. Provided that the /Vs 
depend only on fluxes [i.e. pi = piipj)], because so does JV±p, e \v, A CT , Xg) = \v ± sgn(X a Xs), X a , Xg), 

the physical area, prescription B is uniquely determined 
by the requirement that, for all directions i, the expo- 
nents \XiCi of the holonomy elements A/jr^ (Cj) have a fixed 
constant Poisson bracket with the variable 



(3.5) 



$e\v,X<r,\s) = (47T7^ 1 v / A)3sgn 
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V\ 
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v± sgn(A rT w), A CT ± 



(3.6) 
,Xs), (3.7) 



up to a sign depending on the orientations of the triad 
components. Here, Ipi = \/Gh is the Planck length. Note 
that v is proportional to the volume of the associated 
Bianchi I universe. Then, at least at the level of the 
Poisson bracket algebra, the introduced requirement can 
be understood as the condition that the holonomies pro- 
duce a constant shift in the volume. 

For any of the prescriptions, and owing to the de- 
pendence of pi on the coefficients of the densitized 
triad, the elementary operator J\p, t generates in fact a 
state-dependent non-linear transformation on the basis 
of states \ne,fi(7,fis) — 8>i|m»)- It i s possible to relabel 
this basis with affine parameters instead of the labels fii, 
so that the transformation generated by N^ is simply a 
shift in the labels. Each prescription requires a different 
affine reparametrization. 



2. Quantum representation in case A 

In scheme A, since pi only depends on p i; the three 
fiducial directions are not mixed, and we can calculate 
a new parameter Vi(/j,i) for each fiducial direction to re- 
name the states in our basis, so that the operator N^ 
generates a constant shift in the new label «j, as was 
done in Refs. 0,0] ■ The action of the basic operators on 
the relabeled states \vg,v a ,vs) is given by [3. [ll| 



Jy±fMjv,X a) Xs) = 



p CT |w,A CT ,A«5) = (47T7^ 1 VA)3sgn(A (T )A^|'u,A (T ,A 5 ). 

(3.8) 

The actions of M±a s and p$ can be obtained from Eqs. 
(|3"77|) and (j3~gj) by interchanging X a and X s [H- When 
we construct the Hamiltonian constraint in Sec. IIV1 we 
will see that states with v = — and a fortiori with 
vanishing Ao- or Xs — are removed from our kinematical 
Hilbert space, and therefore the above representation is 
well defined. Unlike in case A, note that all directions 
are mixed now, and the operators Np, 
do not commute. 



and Nfy (i ^ j) 



B. The inhomogeneous sector 

For this sector we employ a Fock quantization, pro- 
moting the variables a m and a* m to annihilation and 
creation operators a m and a\ nl respectively, such that 
[a m , atj = 5 m m- As before, we assume that these oper- 
ators are the identity acting on the homogeneous sector. 
We call S the vector space whose elements are finite lin- 
ear combinations of n-particle states, 



|n) := \...,n-2,n-i,ni,n 2 ,...), 



(3.9) 



Pi\Vi) 
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(671-7^1 



A)3 S ga(vi)\vi\3\vi). 
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with ^2 m n m < 00, n m € N being the occupation number 
(or number of particles) of the 771-th mode. Then, the 
symmetric Fock space T is the completion of the space 
S with respect to the Fock inner product (n'|n) = 5 n ' n . 
Therefore the n-particle states provide an orthonormal 
basis of the Fock space T. 

In the totally deparametrized system, this is the unique 
Fock quantization in which the field dynamics is unitarily 



implemented and that also provides a natural unitary im- 
plementation of the gauge group of S 1 translations J19l |. 



C. Quantum representation of the S 1 symmetry 

Since the generator of translations in the circle, given 
in Eq. (I2.ip . only affects the inhomogeneities, it is the 
same for the Gowdy model in both schemes A and B. Em- 
ploying the above Fock representation and taking normal 
ordering we obtain its quantum counterpart 



Ce =2_j mX ™> X m = &l 



dm a_ /ty ,CL 



m u- m tfc_ m u_ m . 



(3.10) 



m>0 



The n-particle states annihilated by this operator are 
those which satisfy the condition 



y^ mX m = 0, X, 



m ' b rn 



(3.11) 



m>0 



They form a dense set of a proper subspace of the Fock 
space T that we denote by JF V . It is (unitarily equivalent 
to) the physical Hilbert space of Ref. [18[ . 



IV. THE HAMILTONIAN CONSTRAINT 
OPERATOR 

Physical states must also be annihilated by the op- 
erator that represents the (non-densitized) Hamiltonian 
constraint [25[ . We will carry out a process of densitiza- 
tion that will allow us to give an equivalent (and more 
convenient) description in which physical states will be 
annihilated by the densitized version of the Hamiltonian 
constraint. Actually, this procedure is also adopted in 
the quantizations of Refs. [a, El- 



A. Densitization of the Hamiltonian constraint 

We define the subspace of zero homogeneous volume 
states as the kernel of the homogeneous volume opera- 
tor V = <E>i \/\Pi\ which represents the physical volume of 
the (coordinate cell in the) Bianchi I spacetime associated 
with the Gowdy cosmology. Let then Cq = Cbi + C{ be 
the operator that represents the (non-densitized) Hamil- 
tonian constraint for the Gowdy model, where Cbi de- 
notes the Hamiltonian constraint for the Bianchi I model 
and C% is the term that involves the inhomogeneities. 

One first constructs the operator Cbi following LQG pro- 
cedures. When symmetrizing it, one can always adopt 
a suitable factor ordering such that Cbi annihilates the 
subspace of states with zero homogeneous volume and 
leaves its orthogonal complement invariant. On the other 
hand, one can construct the other operator Cj so that 
it inherits the same properties. Therefore the subspace 



of states with vanishing homogeneous volume decouples 
under the action of the full Hamiltonian constraint Cq , 
and one can then remove it from the kinematical Hilbert 
space and restrict the study to its complement. This 
complement will be denoted by %Kin <8> J 7 , where T-Lkiu is 
the completion of 



Cyl s = spanjlue,^,^); v e v a v s ^ 0} 
= span{|u, Ao-, \$); v\ a \s ^ 0}. 



(4.1) 



In the last formula, we have made explicit that only non- 
zero values of A CT and Xg are allowed, even if this is im- 
plicit in the non-vanishing of v = 2Xe\ (7 Xs. 

In general, non-trivial physical states, which are an- 
nihilated by Cg , are not normalizable in Hkih J~. In 
principle, they should belong to a larger space, typically 
the algebraic dual of a suitable dense subspace of this 
kinematical Hilbert space, e.g. the tensor product of 
Cyl s and a suitable dense subspace of T. We will denote 
these states by (ip\. Actually, they can be transformed 
into other states {ip\ by the map 



M = (i>\ 



(4.2) 



which is a bijection in the considered kind of algebraic 

dual spaces. Here, the operator [1/V] represents the in- 
verse of the homogeneous volume and is well defined in 
Cylg (8> T. In Appendix [Bl we provide its explicit form for 
each of the schemes A and B. In both cases, the resulting 
operator is diagonal in our basis of states, annihilates 
the zero homogeneous volume states and is bounded. 
Therefore it can be extended uniquely to the kinemat- 
ical Hilbert space. Moreover, its inverse [1/V] -1 is also 
well defined via the spectral theorem once we have re- 
stricted the discussion to the kinematical Hilbert space 

T-l-Kin^J 7 , because the discrete spectrum of [1/V] in this 
space does not contain the zero anymore. 

The transformed physical states (i/>| are then annihi- 
lated by the (adjoint of the) symmetric operator 



Cr 



Cc 







BI 



c f 



(4.3) 



It is worth noting that the relation between the Hamil- 
tonian constraint of the Gowdy model and its densitized 
version does not involve the volume of the Gowdy space- 
time, but the volume of the associated Bianchi I space- 
time. Hence, the above operator Cq is in fact the quan- 
tum counterpart of the constraint Cq given in Eq. (12.31) . 



B. Quantum representation of the densitized 
Hamiltonian constraint 



For both schemes, we start with the (non-densitized) 
Hamiltonian constraint of the Bianchi I model Cbi • We 



take advantage of the freedom in the factor ordering to 
get a representation as convenient as possible. More 
specifically, as mentioned above we adopt a symmet- 
ric factor ordering which has the two following features: 

i) the same powers of \pi | appear on the left and right of 
every term, in a way such that Cbi decouples the zero 
homogeneous volume states; ii) the factors of the type 

sin(^Ci)sgn(pi) [where sin(^Cj) = i(A/l 2 ^ - 7V r 2/ij .)/2], 
are symmetrized in the form 



sin(^ i c l )sgn(p l ) + sgn(pi)sm(M»Ci) 



(4.4) 



As a consequence, our operator Cbi also decouples states 
with different orientations of the densitized triad com- 
ponents. Both properties have relevant consequences, 
as was discussed in detail in Ref. [5J, where the flat 
Friedmann-Robertson- Walker model coupled to a scalar 
field was quantized adopting the same procedure. As we 
will see later on, the procedure leads to simple superselec- 
tion sectors with neat physical properties. In particular, 
this fact allowed us to solve explicitly the Hamiltonian 
constraint and determine the physical Hilbert space of 
the hybrid Gowdy model for scheme A in Refs. |20ll2l|. 
Therefore, it seems most reasonable to apply the same 
kind of symmetrization in case B as well. 

By densitizing Cbi we obtain then the densitized 
Bianchi I term Cbi- The contribution of the inhomo- 
geneities is contained in the term Q , which is constructed 
by promoting the second line of Eq. (I2.3P to a symmet- 
ric operator. In particular, the comparison between Cbi 
and its classical counterpart [first line of Eq. (I2.3[) ] gives 
us a natural quantum prescription to represent the term 
(c<jp a + cgps) 2 ■ In addition, we also know how to repre- 
sent the term l/|j?g|, as explained in Appendix iBl These 
terms, acting on the homogeneous sector, have a different 
representation in schemes A and B, and we will analyze 
them for each case separately. But let us deal first with 
the terms that have a non-trivial action in the inhomo- 
geneous sector, which are the same for both cases. 



1. Inhomogeneities 

Choosing normal ordering, the quantum counterparts 
of the free Hamiltonian and the interaction term are 

oo 

#0 = X! TO ^ m > N m = a* m a m + aL m a- m , (4.5) 

m>0 



^nt = Yl ~ ~' *™ = & rna- m + O^ol^. (4.6) 



ro>0 



The inhomogeneous sector of the kinematical Hilbert 
space, i.e. the Fock space J 7 , can be written as a di- 
rect sum of dynamically invariant Fock subspaces. In- 
deed, the operator Y m , which is the only operator in the 
Hamiltonian that does not act diagonally on the basis of 



states |n) of J 7 , annihilates and creates pairs of particles 
in modes with the same wavenumber (i.e., the modes with 
wave vectors to and — m). Therefore the quantities X m , 
defined in Eq. (I3.11[) , are conserved under the action of 
the Hamiltonian constraint Cq. Hence, it is convenient 
to relabel the basis of n-particle states with the quantum 
numbers X m , for all positive integers to, together e.g. 
with the eigenvalues N m — n m + n^ m of the operators 
N m , defined in Eq. (|4.5[) . That is, we rewrite the states 
in our basis as 



|A 1 ,A 2 ,...;7V 1 ,A 2 , 



\x-,m). 



(4.7) 



Here, the numbers X m can take any integer value, 
whereas N m <G {|A m | + 2k. k e N}. 

On these states, the action of the relevant operators is 



X m \X;Vl)=X m \X;Vl), 
N m \X;% = N m \X;Vl}, 



(4.8) 
(4.9) 



Y m \X;Vl} = ^™ n X ™ \X;...,N m -2,.. 



y/(N m + 2)2 - XI 



|X;...,A m + 2,...), (4.10) 



and thus the sequence X = {-Xi, A 2 , ...} is not affected, 
as we pointed out above. 

In addition, we will denote by S% the subspace of iS 
spanned by the n-particle states with fixed sequence X, 
and by !F% the respective completion. Then 



T 



^x^x- 



(4.11) 



In practice, as far as the action of the Hamiltonian con- 
straint operator is concerned, we can restrict the study 
of the inhomogeneous sector to any specific subspace Tx- 
In what follows, to simplify the notation, we will denote 
the n-particle states by |9T) and obviate their dependence 
in the fixed sequence X. 

Obviously the operator Hq with domain S% is well de- 
fined in Tn, because it acts diagonally on the n-particle 
states and maps S% into itself. On the other hand, the in- 
teraction term H^ nt creates infinite pairs of particles, and 
thus it does not leave invariant the domain 5^, which 
only contains states with a finite number of them. From 
Eqs. P~oT) . (|3T9l . and (|4~TU)) . we have 

^ \m\ 

X m + 2N m Tj—y 



wL\m 2 



N 2 



ro>0 



2m 2 



(4.12) 



m>0 



Since in the n-particle states only a finite set of occu- 
pation numbers differ from zero, among the above three 
sums only the last one involves an infinite number of non- 
vanishing terms. Actually this sum converges, and hence 
the norm of H^ nt \ 91) is finite. In conclusion, iJ^J^t) € J~x 
for all \y\) € Sx, and therefore H^ nt with domain S% is 
also well defined. 



2. Harniltonian constraint in scheme A 

Following the procedure sketched in the beginning of 
this subsection, for the densitized Harniltonian constraint 
and in scheme A we get an operator of the form (see Ref. 
2l| for the details of the construction) : 



c A - 2 

7 


e 9 e s + e e e CT + e a e s 


+ *,{ (§ - 


7 2 


1 

\-V\M- 


hL + 



327T 2 Qi^ 



(4.13) 



where Hq and H^ nt are given in Eqs 



(1331) and ([HD , 
respectively, \pg\ is constructed from Eq. (I3.4[) . and 

[l/-\/|pe|] is defined in Eq. (jB2[) . The symmetric operator 
Qi represents the classical quantity CiPi, and its action on 
the basis of states \vi) of the homogeneous sector, as de- 
termined by the quantization procedure explained above, 



has the form [11 1 



&i\vi) = -in^[U(vi)\vi + 2) - f-(vi)\vi - 2)]. 

(4.14) 

Here, f±(vi) are two positive functions which satisfy that 
f—(vi) = f+(vi— 2) and, remarkably, that f+(vi) vanishes 
in the whole interval Vi £ [—2,0). Their explicit expres- 
sions are provided in Appendix [Cl 

It is worth noticing that, in scheme A, the homoge- 
neous sector is completely factorized in three indepen- 
dent directional subsectors. 



3. Harniltonian constraint in scheme B 

Motivated by our previous analysis, carried out for the 
Bianchi I model LLpl as well as for the hybrid Gowdy 
model in scheme A [20, [2l[ , we follow in case B the very 
same densitization and symmetrization procedure. 

Let us focus first on the Bianchi I term. Taking into 
account Eq. (I3.1[) for prescription B, one can see that the 
densitized term Cgj corresponds to the symmetrization of 



^ F 2 sgn(p J )sin(^c : ,)sgn(p fc )sin(^ fc Cfe). (4.15) 



7 2 A 



j<k 



This is precisely the gravitational part of the constraint 
represented in Ref. [12], up to a constant multiplicative 
factor (26(. Nonetheless our symmetrization is different. 
Adopting a symmetrization similar to that introduced in 
previous subsections, we get 



--BI 



C = c< 9 » + cm+cw, 



c& 



F 



VlFjVFk+FkVFjWV 



4 7 2 A 
sin(/i l c l )sgn(pi) + sgn(p l )sin(/l l c l ). 



(4.16) 

(4.17) 

(4.18) 



The action of the operators F, is displayed in Appendix 
ID 11 while the action of V = ®»-v/|pi| is obtained from 



Eqs. ([321) and fl37 

The final result for the densitized Harniltonian con- 
straint in scheme B is 



& = c^ + c^ + c^ 



+ 1 



PI 



\po 



G 



\po 



Hi t +32K 2 \p e \Hl}, (4.19) 



where Hq and H^ nt are again given in Eqs. (|4.5[) and 

(|Q|) . [l/|p e |'] can be found in Eq. (|B3l) , \p g \ is con- 
structed from Eq. p.6p . and we choose the symmetric 
operator 



G 



4 7 2 A 



y/V[F a VF a + F 5 VFs]VV- C w (4.20) 



to be the quantum counterpart of [(c a p a + csps)/j] 2 - 

In contrast with scheme A, the homogeneous sector is 
not factorized anymore in three independent directional 
subsectors, since the operators Fi and Fj (i ^ j) do not 
commute. 



C. Superselection in the homogeneous sector 

1. Superselection in scheme A 

As discussed in Refs. [20, 21], Cq leaves invariant the 
Hilbert subspaces rlf., defined as the Cauchy comple- 
tion (with respect to the discrete inner product) of the 
subspaces 

Cyl| ej =span{|^>; Vi e£±}, (4.21) 

where Cf. denotes the semi-lattice of step two defined by 

L% ={±(e l + 2k); fc e N}, e t £ (0,2}. (4.22) 

Therefore, for each £,, the kinematical Hilbert space 



Hi ® J 7 , with ni = ®iH: 



(4.23) 



provides a superselection sector, and the Harniltonian 
constraint operator is well defined in any of its subspaces 
Hi ® J-je, with dense domain ®»Cylg e . <S> S%. Here, we 
are assuming that physical observables can distinguish 
between different modes and thus they do not superse- 
lect Fx- 

In conclusion, in case A we can restrict the study to a 
kinematical Hilbert space that is separable, and in which 
the quantum numbers representing the homogeneous de- 
grees of freedom are strictly positive, with minimum val- 
ues Si, and distributed on cubic lattices of step 2. 



2. Super-selection in scheme B 



and, if v* > 2, also 



Let us analyze the action of the Hamiltonian constraint 
operator Cq on the homogeneous sector component of 
any state \v*,X*, A£)<g)|9T). The explicit expression of this 
operator is given in Appendix ID 21 First, concerning just 
the homogeneous sector, it is straightforward to see that 
the action of the constraint leaves invariant the subspace 
of positive densitized triad coefficients, given by 



Cylg = span-fJ^A^Aa); v,X a ,X s > 0}. 



(4.24) 



We restrict our discussion to this subspace from now on. 

It is worth commenting that, whereas Ref. [12| adopts 
the same symmetrization for the powers of V that we 
have proposed (see e.g. Refs. [20, l21|). so that the 
states with vanishing (homogeneous) volume are indeed 
decoupled, the symmetrization chosen for the signs of the 
triad components differs from ours, and therefore states 
with different orientations of those densitized triad com- 
ponents are not decoupled. Moreover, although the dis- 
cussion in Ref. 12;] is also restricted to the space Cylg , 
that restriction is incorporated there on the basis of the 
symmetry under parity, but the (gravitational part of 
the) Hamiltonian constraint operator does not leave in- 
variant this domain, since it mixes states with different 
orientations. In other words, it is only on the subspace 
of parity symmetric (or antisymmetric) states that the 
restriction of the Hamiltonian constraint to the sector of 
positive orientations can be defined, because the symme- 
try under parity allows one to identify states with nega- 
tive orientations with a counterpart in this sector. 

The Hamiltonian constraint produces shifts on the 
variable v* which are equal to 4 or —4. Nonetheless, 
in the considered sector, a shift by —4 is possible only if 
v* > 4. Then, if we define 



C e = {e + 4k; keN}, £G (0,4], 



(4.25) 



the action of Cq never mixes states with values of v in 
different semi-lattices of this kind. Furthermore, act- 
ing on a state \v*,X*,X$) <g> |9t), this operator produces 
new states with the following quantum numbers v and 
(A a , Xb), where the two identifications (a, b) — (a, S) and 
(a, b) = (S, <j) are allowed: 



K>^prK), K^)- 



v*-2 y* 

<A„, — ■ —At, 



V* 



v* -2 



b > 



• V — V 



X,, . 



2 
4 

v* + 2 



K ) ^O! 



-K 



u^+2 t^+4 A / ^+4 



We see that the effect caused on the A-labels does not 
depend on the reference quantum numbers A* and AJ, 
but only on the value of v* = e + 4k*. This depen- 
dence is through fractional factors whose denominator 
is two or four units bigger or smaller than the numera- 
tor. Therefore, it is possible to see that, starting with 
\v*, A£, A£) <8> \yi) and restricting the consideration to the 
given value v* of the label v, the iterative action of the 
constraint operator leads only to states whose quantum 
numbers A a are of the form A a = uJ e X* (a = er, <5), with 
uj e belonging to the set 



m 




n 



e + 2m 

e + 2n 



(4.26) 



where fc™ € N, and z e Z if e > 2, while z = when 
e < 2. The discrete set W e is countably infinite and turns 
out to be dense in the positive real line. The proof of this 
last statement can be found in Appendix ID 31 Therefore, 
whereas the variable v has support on simple semilattices 
of constant step, the variables A a take values in much 
more complicated sets. Nonetheless, they are also super- 
selected in separable sectors, a fact which had not been 
realized in previous literature. As a particular case, we 
can see that if e = 2 and A* is a fraction, then A a can 
take any value in the set of positive rational numbers. 

In conclusion, Cq leaves invariant the Hilbert sub- 
spaces T-Ls,X" ,\* defined as the Cauchy completion with 
respect to the discrete inner product of 



A* V *-\A X* V *~ 2 X*\ 






c yh,e,\t,\* 5 =span{|«, X a ,X s );v e C e ,X a = u e X* a , 

w e € W e ,A*eM+}. (4.27) 

As a consecuence, the Hilbert subspaces 



A„ , 



" A 6 ' ^oi 1 — K I 



~~^ — a oj ~z t K ) i (K'K) 



u 



e,A* ,At 



•T 



(4.28) 



provide superselection sectors. Moreover, the Hamilto- 
nian constraint operator Cq has well-defined restrictions 
on any of the subspaces H e ,\* ,x* <8> F%, with correspond- 
ing dense domain given by Cyl s E A * A * ® S%. 



V. PHYSICAL SECTOR IN CASE A 
A. Imposition of the Hamiltonian constraint 

In Ref . j2ll] we already determined the solutions (ip\ of 

the Hamiltonian constraint Cq, which were given by an 
expansion of the form 

(ip\ = V] / du a duj s {v e \ ® {u a \® (us\® {il)u a ,u s {ve)\. 

(5.1) 

Here, \u a ) (a — a, S) denotes the generalized eigenstates 
(with real generalized eigenvalue u> a ) of O a . These oper- 
ators are constants of motion (they commute with Cq), 
and therefore the states \tu a ) are stable under the ac- 
tion of Cq. The concrete expression of these states can 
be found in Ref. [ll| . On the other hand, remarkably, 
the projections (VW-^,s( e + 2/c)| of the solution on the 
sections with constant value of vg = eg + Ik for every 
k 6 N + , are all formally determined by the initial data 
{i>u a ,ui s (Ee)\ through the action of a complicated opera- 
tor, which involves the iterative action of Hq and -ff int . 
Again, the explicit expression is provided in Ref. [2l| . 

In addition, the physical solutions must be annihilated 
by the generator of translations in S 1 . This S 1 symmetry 
is preserved by the dynamics, since the quantities X m are 
constants of motion; therefore we can ignore it at this 
stage and impose it after dealing with the Hamiltonian 
constraint. 



is an extremely difficult task, given the complexity of our 
model, and no satisfactory choice is actually at hand. In 
this sense, the kinematical structure of our quantization 
is not well adapted to the physical one, and their relation 
is not straightforward a priori. 

Nonetheless, we can still complete the quantization 
program: we just need to make sense of the solutions 
and provide them with a Hilbert structure. As we have 
already said, the solutions to the densitized Hamiltonian 
constraint are completely determined, at least formally, 
by a single piece of initial data (Vw,wa(£fl)|> and then we 
can identify these solutions with the corresponding data. 
The determination of a complete set of real classical ob- 
servables acting on these initial data, together with the 
condition that they be represented as self-adjoint opera- 
tors, determines a unique inner product [28j | that charac- 
terizes the Hilbert structure. Moreover, taking into ac- 
count the additional S* 1 symmetry, implemented by con- 
dition (I3.11[) , we conclude [2CJ, |2l| that the initial data 
Ww,a;j(£0)| must belong to the Hilbert space 



L z 



,du a dus) ® T v 



(5.2) 



that we identify as the physical Hilbert space. 

Alternatively, we can argue that this is the physical 
Hilbert space following a different line of reasoning, based 
on the idea that we can regularize the theory by means of 
a cutoff for the wavenumber m > and analyze the limit 
of arbitrarily large values of the cutoff. Note that we are 
allowed to carry out this further reduction of the phase 
space because, in the Gowdy model, the field modes with 
different wavenumbers are not mixed dynamically. 



B. Physical Hilbert space 

The resulting form of the general solutions (-01 of the 
Hamiltonian constraint is only formal, in the sense that 
they do not belong to the dual of the domain of definition 
of Cq. Namely, some of the coefficients (Vw.^f (^fl)]^) 
of the solutions diverge, for example when vg > eg + 4. 
Indeed, it is not difficult to see that the space Sx is not in 
the domain of the operator HqH^, because this contains 
the term J2 m ^mXm- The action of this term on a generic 
state of Sx does not lead to a normalizable state of the 
Fock space Tx, not even in the (conventional) generalized 
sense, since, in particular, it involves the creation of pairs 
in an infinite number of modes and then a sum over the 
number of all of those particles. 

This problem can be traced back to the choice of do- 
main for the Hamiltonian constraint, since the present 
one is not invariant under the action of Cq. Indeed, it is 
worth pointing out that the determination of an invari- 
ant domain (and such that the Hamiltonian constraint 
be an essentially self-adjoint operator) would allow one 
to resort to group averaging techniques [27| in order to 
construct the space of physical solutions. However, the 
selection of an alternative domain that remains invariant 



1. Regularized model 

Let Sm be the subspace of S spanned by the n- 
particle states which satisfy the condition N m = for all 
m > M > 0, and let S M be its orthogonal complement. 
In addition, let Pm and P^ be the projectors on the sub- 
spaces Sm and S M , respectively. We are now interested 
in finding the physical Hilbert space of the truncation 
of the Gowdy model in which all the field modes with 
wavenumber bigger than the cutoff M vanish. Obviously, 
this truncated system is governed by the constraints 



(^,u> s (v e )\P M = 0, 
(ipu a ,u s (ve)\PM Cg P m = 0. 



(5.3) 
(5.4) 
(5.5) 



Simultaneous solutions to the first two equations have 
the same form as the solutions to the Hamiltonian con- 
straint in the full Gowdy model, but containing exclu- 
sively wavenumbers m equal or smaller than M, so that 
they now possess only a finite number of terms. As a con- 
sequence, the divergences caused by the infinite produc- 
tion of pairs of particles disappear, and the coefficients 
(.ipw„us(ve)\yi) of the solutions to Eqs. (|S"3|) and ([Q]) 
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are now finite for all vg. Therefore, the solutions (for 
fixed vg, u> ai and cug) live in the dual S^ of the space 
spanned by the n-particle states with the cutoff. Simi- 
larly to what we did in the previous subsection, we can 
now endow the data at vg — eg with a Hilbert structure. 
We construct the same complete set of observables [21| , 
with the difference that now we do not have an infinite 
number of them in the inhomogeneous sector, but only 
AM (corresponding to two degrees of freedom on phase 
space for each of the 2M wave vectors). The resulting 
observables act on the solutions (Vw,w 5 ( £ 0)l and are self- 
adjoint in L 2 (R 2 ,duj a doJs) <8> J~m, where Tm is the Fock 
space obtained by completing Sm- Taking into account 
the remaining constraint (|5.5|) . which implements the S 1 
symmetry, the physically admissible states have finally 
the Hilbert structure 



,duj a dtus) <g> (.7>)m, 



(5.6) 



where (-F p )m is the subspace of Tm spanned by the 



M 



0. 



particle states that verify the condition X)m>o rn -^- 

As the cutoff M increases, we get closer to the non- 
truncated theory for the Gowdy model. This indicates 
that, in the limit M -> oo, we should recover the full 
Fock space: (J- p )m — > J~ p . This supports our previous 
statement that the physical Hilbert space of the non- 
truncated Gowdy model is in fact the space (|5.2[) . 

The introduction of the cutoff in the wavenumber also 
makes manageable the numerical study of the effect that 
the inhomogeneities have on the dynamical behavior of 
the system and, in particular, of the changes that oc- 
cur in the dynamics of the Bianchi I background when 
the inhomogeneities grow. Actually, this kind of analysis 
has already been carried out at the effective level and 
is included in the study of Ref. [29|], where the classi- 
cal Hamiltonian has been modified in an effective way 
in order to take into account the corrections that arise 
from the quantum theory — although the results obtained 
there are extended to the full model without the cutoff. 



VI. PHYSICAL SECTOR IN CASE B 



A. Imposition of the Hamiltonian constraint 



Namely, 

®(V(«,w e Aj,w e AJ)|. (6.1) 

Inserting this expansion into the (dual of the) con- 
straint equation, projecting on the homogeneous sector, 
and taking into account the action of Cq given in Ap- 



pendix [1721 we g e t that the projections 

(i(){v,\ a ,\s)\ = (ip(v,uj s Xl,iD e X* s )\ 

satisfy a series of relations that can be interpreted as 
difference equations in v. Introducing the projections of 
(■01 on the combinations of states defined in Eqs. (|D5|) - 
(|TJ8|) of Appendix ID 21 which we call 



(V>±(« ± 4, A CT , X s )\ = (i;\v ± 4, A CT , X s )± 

and in a similar way for the rest of projections, the rela- 
tion obtained is 

(ip+(v + 4, A CT , X s )\ - /3b 2 e (v, A CT , X s )b 2 g (v + 4, A CT , X s ) 

x^0A> + 4,A ff ,A,)|#L 
1 32v 2 



-pXlX 2 x + {v) { ^ X ^ X5) ^ 



■ r ' ; (W W («, A., \s)\ + ^4(^o + (v, A., X S )\ 



x+(v) 

X-(v) 

x+(v) 



x+(v) 
{i/>-(v - 4, X a ,X s )\ + (3b 2 g (v, X a , X s ) 



■ { b 2 6 (v - 4, X a , X 5 f 4 X ~^\ W_(v - 4, K, X 5 )\ 
1 v x+(v) 



6g(u, Ao-jAj) 

x o( v ) , i * 



' ,T,y) (^-(«,A CT ,A 5 )| 



x+(v) 



x+(v) 
(ip' 0+ (v,X a ,X s )\ 



Hi 



(6.2) 



Here, b$, x±, and x are the functions defined in Appen- 
dices E2] and EU 



Unlike in case A, the homogeneous sector of the kine- 
matical Hilbert space, T-L e ,\* ,\* > is not factorized in three 
independent directional subsectors, and none of the op- 
erators vVFiW appearing in the expression of the con- 
straint Cq (i.e., the counterpart of 0; in case A) is now 
a constant of motion [3Pj. Therefore, in this case we 
cannot simplify the action of the Hamiltonian constraint 
operator by diagonalizing it in the directional subsectors 
labeled by a and 6, as we did in scheme A. 

Then, in order to solve the Hamiltonian constraint 
represented by Cq, we expand the solutions (4>\ using 
the basis of states \v, X a , Xs) for the homogeneous sector. 



B. Analysis of the solutions 

We want to prove now that, as in case A, the solution 
is completely determined by the data on the first sec- 
tion v = £ (at least formally). Specifically, we want to 
show that, given a set of initial data (i/)(e, U) E X*, w E A|)| 
(belonging to S% for all w e , Co e £ W E ), it is possible to de- 
termine each term (ip(v, X ai Xs)\ of the solution, for every 
v > e in L e . 

The presence of the interaction term on the left hand 
side of Eq. (16.21) complicates a direct proof of the above 
statement. However, it is possible to attain the result 
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by means of an asymptotic analysis of the solutions. Re- 
markably, our theory involves a dimensionless parame- 
ter, /3, introduced in Eq. (|D4p . and recurring to it we 
can naturally adopt an asymptotic approach without the 
need to introduce any external parameter by hand. Note 
that, since the area gap A is proportional to 7^, /3 is 
proportional to the inverse of the Immirzi parameter 7. 
Thus, in the limit j3 — > 0, we expand the solutions in 
asymptotic series of the form: 



(V(£ + 4*, A CT , Xe)\ = J2 0"~* (X £ + 4fc > x -> x s)\ , 
Vfe € N+. (6.3) 



Note that the linear combinations introduced in Eq. 
(|6.2p . like e.g. (ip±(v ± 4, A CT , Aa)|, adopt then similar 
expansions, and we will denote their terms using an ob- 
vious notation, for instance ("">p±(v ± 4, A CT , A<s)j. Sub- 
stituting the expansion (J6.3I) in the constraint (|6.2p . and 
considering powers of /3 order by order, we obtain an ex- 
pression for every term ( n i(j + (v + 4, A CT , A^)| (for generic 
w) provided that the data at v and u — 4 are known. The 
explicit result is the following (to simplify the notation, 
we obviate the dependence of the states and of the func- 
tion be on the A's): 



Leading term: 



(>+(t> + 4)| = - 



32v 2 



\l\lx + {v) 



(V(t/)|flg, (6- 4 ) 



• first order correction: 
32 V 2 



(V+(« + 4)|=- 



^^x+(«) 



(V(v)|fl§ 



H^4(Vo-(«)| + ^S(Vo + W] 



x+(v) 



ar+(u) 



+ $(t,)$(t; + 4)^(y + (t, + 4)|fl£ t , 



(6.5) 



• second order correction: 
32w 2 



V+(« + 4) 



A 2 A 2 x + ( v ) 



(V(«)|^ 



+^(vo-(.)l+^t 1 

«+(«) x+(u) 



^o+( u )| 



^r(V-(--4)|+fo^){-6 2 («) 



.2, , „^ + 4 n,/.' 



+ &S(v + 4)- 



^;c«+4) 



H, 



(6.6) 



• n-th order correction (n > 3) 
32^ 2 



( n ^+(v + A) 



AJAKM ( "* ( " )|ff ° 



a; + (w) v [ 

x+(u) ' x + (v) 

+ ^( t , + 4 )l±i(»-V + ( t , + 4)| 



+ 6 2 ( U -4)^-i("" 3 ^( v -4)||>//;; !i . ((;.',) 



The above expressions simplify considerably in the 
cases v = e and v = e + 4. This is due to the fact 
that, on the one hand, the data on the very initial section 
are not given by asymptotic series, that is, in principle 
(^(e)| = (°ip(s)\ and hence ("V0)l = for all n > 1. 
And, on the other hand, X—(e) = for all e <G (0,4], 
whereas x (e) = if e < 2 [see Eqs. (JD9J) and (|PTT|) ]. 

In view of these equations, we see that, for v > 4, 
the knowledge of the solution on the sections v — 4 and 
u, together with the terms of asymptotic order n — 1 
on the section i> + 4, determine on this last section the 
terms of the solution at order n via the following linear 
combinations 



("V+(^ + 4,A ff ,A 5 )| = "V [v + i,X, 



v + A 
v + 2 



X S 



+ [ n ip[ v + 4,\ a ,^^X s 



. .. . . v + 4. v + 2 
+ [ n ip [v + 4, — — A ff , X s 

v + 2 v 



+ [ n ip[ v + 4, ^-^A CT ,A 5 



+ 'V U + 4 



-y + 2 

v + 2 , v + 4 



-Xn 



v + 2 



Xs 



+ (™V( v + 4, ^— A CT ,A 5 



(6.8) 

On the other hand, for < v < 4, the data on the section 
v — 4 are spurious, and their knowledge its not required 
to determine the terms of the solution for v + 4. 

Actually, a similar structure appears also in the so- 
lutions to the Hamiltonian constraint for the Bianchi I 
model, inasmuch as the solution on the section v + 4 is 
determined in terms of the same kind of linear combina- 
tions. Furthermore, it has been recently shown [31] that, 
for all v > 0, the set of linear combinations 

{(<ip+(v + 4,uj s X*,u>eX* s )\; We,w e e>V £ ,} 
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determines the set of individual terms 

{(^(» + 4,w e Aj,w e Aj)|; u e ,(D e €W e } 

through a one-to-one map. This result can be applied 
here as well as for Bianchi I cosmology. Therefore, start- 
ing with the initial data, we can obtain, step by step, 
the terms of the solution — up to the desired asymptotic 
order — on all of the consecutive w-sections. 

In conclusion, the initial data {ij){e, A CT , Xg)\ (with A CT 
and Xs taking all possible values in the corresponding su- 
perselection sectors) completely determines the solution, 
as we wanted to show. The solutions, constructed in this 
way, are formal, as happens to be the case in scheme A, 
in the sense that the objects ( n ip(v + 4, A CT , A^)| do not 
belong in general to the dual space of the domain chosen 
to define the Hamiltonian constraint operator, owing to 
the presence of the operator W L in their expressions. 



quantized Bianchi I background in the space of physical 
states, starting with a hybrid quantization in the kine- 
matical setting. Indeed, we have proved that this is the 
case, since the physical Hilbert space obtained in both 
schemes has the structure of the tensor product of the 
physical Hilbert space of the Bianchi I model and a Fock 
space, which turns out to be equivalent to the space ob- 
tained in the standard Fock quantization [18|, [l9[ . This 
result supports the validity of the hybrid quantization, 
because the latter should lead to the standard quantiza- 
tion of the system in the limit in which the effects arising 
from the discreteness of the geometry become negligible. 
Let us remark that the result is non-trivial, inasmuch 
as the hybrid approach is introduced in the kinemati- 
cal arena and the relation between the kinematical and 
physical Hilbert structures cannot be anticipated before 
completing the quantization, even more if one takes into 
account the field-like complexity of the model. 



C. Physical Hilbert space 

In order to provide the solutions with a Hilbert struc- 
ture, we proceed in the very same way as in case A. Once 
we have shown that the set of initial data 

{(ijj(e,uj e Xl,Ld e X* s )\ ; u e ,u e £ W s } 

characterizes the solution, we can identify solutions with 
their corresponding data, and the physical Hilbert space 
with a Hilbert space of such initial data. 

The requirement that a complete set of observables 
acting on these initial data be self-adjoint operators de- 
termines again uniquely the inner product that provides 
the Hilbert structure. Such observables are given, for in- 
stance, by a complete set of observables for the Bianchi 
I model in vacuo and by the observables introduced in 
scheme A for the inhomogeneities. Imposing the remain- 
ing S 1 symmetry on the resulting Hilbert space, we fi- 
nally get the same structure found in case A for the 
physical Hilbert space, namely, the tensor product of the 
Fock subspace T v and the physical Hilbert space of the 
Bianchi I model in vacuo, though now in scheme B: 



^phys — ^phys.BI ® ->p- 



(6.9) 



The explicit form of "Hp hvs BI is analyzed in Ref . [3l| . 



VII. SUMMARY AND CONCLUSION 

A. Recovery of the standard quantum field theory 

As we have commented in the Introduction, one of the 
motivations of this work is to investigate the plausibility 
of the recovery of standard quantum field theory in the 
framework of loop quantization. In particular we wanted 
to show that, in the Gowdy model, one attains a Fock 
description of the inhomogeneities over a polymerically 



B. Resolution of the cosmological singularity 

The classical solutions of the linearly polarized Gowdy 
T 3 model present generically a cosmological singularity. 
In Ref. [15| . e.g., a curvature invariant was explicitly 
calculated and proven to diverge almost everywhere at 
initial time. In terms of the variables that we have em- 
ployed, the cosmological singularity corresponds to van- 
ishing values for the components pi of the densitized 
triad. Actually, as we can see in Eq. (1A1[) . the metric is 
ill defined if any of the p^s is zero. 

In our quantum theory, the polymeric quantization 
performed in the homogeneous sector succeeds in elim- 
inating the singularity. More explicitly, we have been 
able to remove the kernel of all the operators pi and, as a 
consequence, an analog of the classical cosmological sin- 
gularity does not exist any more quantum mechanically. 
This resolution of the singularity is achieved at a kine- 
matical level. Of course, it persists in the physical Hilbert 
space, since physical states do not have projection onto 
the zero eigenspaces of the operators pi. Furthermore, 
they only have support on a sector with fixed orientation 
of the triad components and, then, they do not cross the 
singularity to another branch of the universe correspond- 
ing to a different orientation. 

On the other hand, in addition to this kinematical res- 
olution of the cosmological singularity, it is worth com- 
menting that, at least for scheme A and in the framework 
of the effective description corresponding to the hybrid 
quantization put forward here, the numerical simulations 
performed so far for the Gowdy model show the pres- 
ence of a bounce which replaces the singularity and which 
emerges owing to the quantum geometry corrections to 
Einstein's theory [29( . Similar numerical calculations are 
being developed currently for scheme B [32|, in order to 
validate this stronger result about the resolution of the 
singularity. 

Let us emphasize that the standard (non-polymeric) 
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quantum methods do not succeed in resolving the cosmo- 
logical singularity. On the one hand, in the Fock quan- 
tization of the deparametrized system, which has been 
discussed in the literature |17[ and which has been suc- 
cessfully accomplished till completion in Refs. [18|, [l9(, a 
classical time parameter is present explicitly in the quan- 
tum description, and the curvature invariant calculated 
in Ref. 15] depends on its inverse in such a way that the 
invariant still blows up at initial time. On the other hand, 
if one does not deparametrize the system, following our 
gauge reduction, and quantizes the homogeneous sector 
in a standard (non-polymeric) way, as in the Wheeler-De 
Witt approach, then the zero eigenvalue would be in- 
cluded in the continuous spectrum of the triad operator, 
instead of the discrete spectrum, and there would not 
exist a proper subspace associated with this eigenvalue 
which could be decoupled and removed. 



C. Concluding remarks 

In conclusion, we have rigorously constructed a hy- 
brid quantization of the Gowdy model with three-torus 
topology and linearly polarized gravitational waves. The 
homogeneous sector of the phase space, which coincides 
with the phase space of a Bianchi I model, has been 
polymerically quantized, whereas we have applied a (dis- 
tinguished) Fock quantization to represent the inhomo- 
geneous sector. In the LQC literature, there exist two 
different schemes for the polymeric quantization of the 
Bianchi I universes, denoted in this paper as schemes A 
and B. We had already analyzed the hybrid quantization 
of the Gowdy model adopting scheme A in Refs. |20l.l21|. 
Here, we have revisited that quantization and extended 
our hybrid approach to the alternative case B, in which 
the homogeneous sector has a different representation. 

In both schemes, the quantum Hamiltonian constraint 
has been densitized, in order to deal with a simpler con- 
straint, and then has been promoted to an operator, well 
defined in some dense domain of the kinematical Hilbert 
space. This is truly a non-trivial result, because our sys- 
tem possesses an infinite number of degrees of freedom 
and the two sectors, on which the constraint operator 
acts, are coupled and quantized with entirely different 
methods. 

As we have seen, the kinematical structure over which 
we have defined the theory and the choice of domain for 
our quantum operators do not suffice to make sense of the 
formal solutions to the Hamiltonian constraint. Nonethe- 
less, we have found a procedure to overcome the problem 
and complete the quantization. Indeed, in our hybrid ap- 
proach, the Hamiltonian constraint provides a difference 
equation in an internal discrete parameter (vg in scheme 
A, v in scheme B) which has a strictly positive minimum 
value, and the solutions to the Hamiltonian constraint 
are completely determined by the data provided on the 
initial section of such a discrete parameter. One can say 
that the solutions follow a no-boundary prescription, in 



the sense that they arise in a single section without the 
need to impose any particular boundary condition. This 
behavior has two important consequences. On the one 
hand, this immediately resolves the classical singularity 
in the quantum theory at a kinematical level. On the 
other hand, it allows one to deal with these solutions by 
identifying them with initial data. In this way, we have 
been able to characterize the physical Hilbert space in 
both schemes A and B. Remarkably, as we have pointed 
out, this procedure leads to the recovery of the standard 
quantum field theory for the inhomogeneities. 
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Appendix A: Classical metric 

To derive the form of the classical metric of the lin- 
early polarized Gowdy T 3 model, one can start with its 
expression in the (field) parametrization of Ref. [18[ , ap- 
ply the gauge fixing procedure, and perform a canonical 
transformation from the elementary variables chosen for 
the homogeneous sector in that parametrization to the 
corresponding Ashtekar variables {ci,pi}. A careful cal- 
culation shows that, in our variables, the non- vanishing 
components of the induced three-metric are 
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Besides, owing to the homogeneity of the shift function 
N e , required by the gauge fixing, we can reabsorb it by 
means of the following redefinition of the coordinate 8: 



9+ / dt' N (t')^6, 
Jti 

where ti is any initial time. Then, the spacetime metric 
becomes 

ds 2 = -q ee (^P) N 2 dt 2 + q 0e d6 2 + q aa d<r 2 + q ss dS 2 . 

Here, N is the densitized lapse function, which in our 
gauge is spatially homogenous. The metric of the Bianchi 
I spacetime is the result of ignoring the inhomogeneities 
in Eq. (JXTJ, setting f = = (. 



and similarly for [l/|p CT I 1 / 4 ] and [l/lp^ 1 / 4 ]. Their action 
on our basis of states turns out to be 



\v, X a , X S ) = h f V -\ >* a ^\ \v, A CT , A,), (B4) 



(Att^VA' 



where (a = a, S) 



b e {v,X <7 ,\8)^V2\XaXs\ y/\v + l\-y/\v-l\ 
b a (v,X a ,Xs) = 



V 



^TT\-^/J^T\ . (B5) 



The inverse homogeneous volume operator can then be 
represented as the regularized operator 



Appendix B: Inverse homogeneous volume operator 

Following the procedures of LQG, the classical expres- 
sion of l/(bi| 1_r ) in LQC, for r > 0, is represented by 
the regularized operator (see e.g. Ref. Q) 



\Pi 



(B6) 



II— r 






(Bl) 



where [1//Uj] is the quantum counterpart of the expression 
given in Eq. (|3.ip for each of the considered schemes. 
The choice of r is arbitrary. 



Appendix C: Details of the quantum model for 
scheme A 



For completeness in the presentation, we include here 
some details about the operator 0^ that appears in 
the expression (|4.13j) of the densitized Hamiltonian con- 
straint operator Cq for scheme A. The definition of 0j in 
terms of holonomy and fluxes operators is [IJ, [2l| : 



1. Scheme A 

In case A, and for r — 1/2, one obtains [9|, LLl| 
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Qi = 
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Wi 



\pi 
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(CI) 



Its action on the states \vi) is given in Eq. (J4.14I) . where 
the functions f±(vi) are: 



2. Scheme B 



In case B, choosing r = 1/4 and combining all the 
powers of \pg | , we obtain the expression 
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f±(vi) = g(vi ± 2)s±(vi)g(vi), 



(C2) 



with 



s±(vi) = sgn(«i ± 2) + sgn(^ 

9{vi) 

5 (0) = 0. 
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Appendix D: Details of the quantum model for 
scheme B 



1. The operators Fi 



Expressing each operator Fi in terms of M±jj i , 



straightforward calculation shows that 



n, , x \ .sgn(A CT A«5) -^ 
F \v,\ (r ,Xs) =t 2^ l 



i=+i,-i 



sgn(w) 



+ sgnjz; + 2lsgn(X a \ s )} 



\v + 2lsgn(\ a \ s ), X a , X s ), 

(Dl) 



F a \v, X a ,Xg) = i 



sgn(A 



+ 1 



A„ 



v + 2lsgn(vX <T ),X a + 2l 



v 



sgn{M + 2Zsgn(A CT )} 
,X S ). (D2) 



The action of Fg is similar to that of F a , interchanging 
the roles of A CT and A^. 



2. The operator Cg 

The action of the constraint operator on the states 
\v, X a , Xs) <8> 1 91) of our basis, with v, A CT , and Xs being all 
positive, is the following: 



FB 



C%\v,Xa,X 5 )®\yi) 



(nl^) 2 
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(D3) 
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(D4) 



and we have introduced the following notation: 



\v±4,X a ,X 5 )-_ 
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\v, Xa, Xs)' ± = 2\v, Xa,Xs) + 



v x v±2^ 

V, — -Xa, A,5 



v±2 , w , 
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v v ±2 
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(D8) 



and 



x-(v) = 2^{v - 2)Vu - 4[1 + sgn(w - 4)], (D9) 

cc+(v) =x_(v + 4), (D10) 

»o(«)=2(«-2)t;[l + sgii(v-2)], (Dll) 

x+{v)=Xq{v + 2). (D12) 

3. Support of the anisotropics 

We want to prove that the set W £ is dense in R + . For 
this, we will show that its subset U e , defined as 

( e + Am 1 

W E DU s = l—^;m,neN\, (D13) 

is already dense in the positive real line. 

Let a and b be any two positive real numbers, such 
that b — a > 0. Besides, we define the set 



V s = || +n: ne n| 



(D14) 
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Then there always exists a number s = e/4 + rii e V e 
such that 1 < s(b — a), or equivalently 



1 < sb. 



(D15) 



Let us denote by t = e/4 + ra\ the largest number in V £ 
which is smaller or equal than sa + 1, that is 



sa < — + mi < sa + 1. 
4 



(D16) 



Equations (|D15|) and (|D16|) imply that sa < e/4 + mi < 
s6. These inequalities can be written equivalently as 



a < 



Ami 



•4m 



<6. 



(D17) 



In conclusion, given any two positive numbers a and b 
with a < b, there always exists a number u G U e such 
that a < u < b. As a consequence C/ e , and W £ a fortiori, 
are dense in the positive real axis, as we wanted to prove. 
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